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OS ■ ABSTRACT 

The deprojection of the surface brightness of axisymmetric galaxies is indeterminate 
. unless the galaxy is seen edge-on. In practice, this problem is often circumvented by 

Q ' making ad hoc assumptions about the density distribution. However, one can redis- 

, tribute the density and still project to the same surface brightness. This is similar to 

adding so-called konus densities to the assumed density distribution. 
0\ ' In this paper we investigate the maximum konus density that one can add to 

elliptical galaxies. In particular we focus on the uncertainties in the central densities 
of axisymmetric, elliptical galaxies due to the non-uniqueness of the deprojection. 

For Stackel potentials a sufficient condition for positivity of the phase space dis- 
tribution function exists, which is used as criterion to determine the maximum konus 
density that one can add to a perfect oblate spheroid. For small inclination angles we 
find an uncertainty in the central density of up to a factor two. 

Elliptical galaxies in general have a central density cusp. We therefore also inves- 
tigate the maximum konus densities of cusped ellipticals. For these models we use an 
^sO , approximate criterion, which we have tested on the perfect oblate spheroid models. 

0\ ' For sufficiently small scalelengths, the central density of the maximum konus density 

"*^h \ that can be added to a cusped elliptical is very high. In order to estimate the dynam- 

O ■ ical influence of konus densities on the central region, we calculate the mass fraction 

I 1 Mkon/Mgaj they can add to the center. We show that this mass fraction is at most a 

few percent. We also investigate the dynamical effect of cusped konus densities that 
have p(r) oc r~ a at small radii. We show that konus densities can only be moderately 
cusped (a < 1), and that an increase in cusp slope a, results in a decrease of the mass 
fraction added to the center by the konus density. 

We illustrate all this by the specific example of M32, and show that the uncertainty 
in the central mass due to deprojection is negligible compared to the inferred mass of 
the central black hole (BH) in this galaxy. 
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1 1 INTRODUCTION 

In recent years considerable progress has been made in the 
construction of axisymmetric dynamical models of galaxies 
(e.g., de Zeeuw 1994). Comparison of such models with ob- 
served kinematics has resulted in strong indications for the 
presence of a nuclear black hole (BH) in a number of galaxies 
(e.g., Kormendy & Richstone 1995). In order to construct 
such models one needs to recover the density distribution 
from the observed surface brightness distribution. However, 
Rybicki (1986) showed, by use of the Fourier Slice Theorem, 
that the deprojection of axisymmetric bodies is indetermi- 
nate unless the inclination angle i = 90° (i.e., the galaxy is 
seen edge-on) . In practice this problem is often circumvented 
by making an ad hoc assumption about the density distri- 
bution, e.g., considering only densities stratified on similar 
concentric ellipsoids. However, the Fourier Slice Theorem 



shows that any density distribution whose Fourier Trans- 
form is only non-zero inside a cone with half-opening angle 
90° — i and aligned with the symmetry axis of the Fourier 
Transform of the density distribution (the so called 'cone 
of ignorance'), projects to zero surface brightness. Following 
the nomenclature of Gerhard & Binney (1996, hereafter GB) 
we will call such densities 'konus densities'. Franx (1988) 
gave an example of such a konus density which, when added 
to a triaxial Stackel model, projects to zero surface bright- 
ness. GB investigated the effect of adding or subtracting a 
specific family of konus densities Pk(R, z\io) to a density dis- 
tribution stratified on similar concentric ellipsoids po(R, z). 
Although all densities 

P (R, z) = p (R, z) + f p k (R, z\i ), (1.1) 

project to the same surface brightness if i < io, the under- 
lying density distributions (and therefore the corresponding 
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dynamics) can differ considerably. 

In this paper we investigate to what extent one can add 
or subtract konus densities. In particular we focus on the 
uncertainty of the central density of axisymmetric galaxies 
due to the non-uniqueness of the deprojection. The kine- 
matical evidence for the presence of nuclear BHs arises from 
the finding that the central density is too low in order to 
explain the observed kinematics. This density distribution 
is generally derived from the deprojection of the observed 
surface brightness, which requires assumptions on the in- 
clination angle and the density distribution itself (i.e., the 
density is stratified on similar concentric ellipsoids). We ap- 
ply our analysis to M32, in which the presence of a BH 
has been claimed. The mass of the BH that is inferred is 
of the order of 20% of the total mass inside 1". We inves- 
tigate whether adding konus densities can redistribute the 
central mass such that we can elude the requirement of a 
BH in M32. We show that, although the maximum konus 
density that one can add can have very high central density 
(the konus density can even be cusped), the actual mass it 
adds to the nucleus (and therefore its dynamical influence) 
is negligible compared to the inferred BH. 

This paper is organized as follows. In Section 2, we de- 
scribe the criterion we use to determine the maximum konus 
density one can add to a certain galaxy model. Section 3 
deals with the specific konus densities we use for our in- 
vestigation. In Section 4 we discuss the uncertainties of the 
central density due to the freedom one has in redistributing 
the density distribution. In Section 5 an approximate crite- 
rion is discussed, which is applied to cusped galaxy models 
in Section 6. We discuss the application of all this to the 
specific example of M32 in Section 7. Finally, in Section 8, 
we summarize and discuss our results. 



2 2 THE CRITERION 

Since konus densities have regions of both positive and neg- 
ative density there is a /max such that for / > /max their 
will be regions where the total density p = p ga i + / pkon 
is negative. Therefore, / max expresses the maximum konus 
density under the physical criterion that p > 0. However, 
this criterion ignores whether or not such a galaxy can ac- 
tually be built from its orbit building blocks. In general, 
for / = / max there will be points (_Ro, \ z o\) in the meridional 
plane with zero density. It is easily seen that such models can 
not be physical. Orbits in an axisymmetric potential densely 
fill a limited area inside their zero-velocity curve. For thin 
tube orbits, this area shrinks to a one-dimensional curve. 
Let those curves be described by ((R, z) and parameterized 
by (f>, such that <f> — a in the equatorial plane, and 4> = b > a 
along the symmetry axis R = 0. It is easily seen that it is 
impossible to construct a model that has only one point in 
the first quadrant of the meridional plane (Ro, zo) = (Co, 4>o) 
with zero density, and p > everywhere else. This is illus- 
trated graphically in Figure 1, where three orbits are plotted 
in the meridional (R, z)-plane. The thick, solid curve is the 
curve with C = Co, an d the solid dot is the point (Co, <fio)- As 
can be seen, all three plotted orbits are just allowed since 
they 'avoid' the solid dot, and therefore do not contribute 
density to (Co,<Ao)- Each orbit that contributes density to a 
point (Co,0')i with 4>' > 00 also contributes density to all 
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Figure 1. The grid are curves of constant C and <f>, which are 
the curves in between which orbits in an axisymmetric potential 
are enclosed. The dashed regions are three such orbits. The thick, 
solid curve is the curve with C = Co an d the solid dot is the point 
(Co, <t>o) ( see text). As can be seen, having zero density in (Co, </>o) 
results in zero density along the entire curve C = Co with <j> > </>o . 

points (Co,0 < 4>')- As a consequence, the only way to dis- 
tribute orbits so that p(Co, 4>a) = is by having zero density 
along the entire curve with (Co, 4> > <t>o)- 

In general this criterion is difficult to make quantitative. 
For Stackel potentials, however, the third integral is known 
analytically and (C, 4>) are the prolate spheroidal coordinates 
(A, v) (see Appendix A). Bishop (1987) formulated the prin- 
ciple outlined above by proving that a sufficient condition 
for the distribution function f(E,L z ,Iz) of a Stackel model 
to be non-negative and hence to correspond to a proper equi- 
librium model is that the density should monotonically de- 
crease along lines of constant A; i.e., 

< 0. (2.1) 

In principle, one can still build physical models (i.e., 
f(E,L z ,Li) > 0) that do not obey this criterion (2.1), since 
it ignores the fact that orbits do not distribute density uni- 
formly over the area in the (A, v) plane occupied by that or- 
bit. The criterion, although sufficient is therefore not strict, 
and will lead to an underestimate of the maximum konus 
density that can be added to a Stackel potential. 

As a specific example of a Stackel model we will focus 
on the perfect oblate spheroids (see Appendix A), and use 
criterion (2.1) to investigate the maximum konus density one 
can add to such a density distribution. We note that if one 
adds a konus density to a Stackel potential it no longer is a 
Stackel potential, and therefore the orbits will no longer be 
confined by curves of constant A and v. However, as a first 
order approximation, we can still use the criterion (2.1). As 
we will show, the results presented in this paper are only 
mildly dependent on the actual curve along which we de- 
mand the gradient of the density to be negative (see Section 
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5). 



3 3 THE KONUS DENSITY 

In the following we will concentrate on ellipticals that have a 
luminous density that decays as r~ 4 . The surface brightness 
of such systems falls of as RT 3 at large radii, which is in good 
agreement with observations. In order for the total density, 
P = Pgai + / Pkon, of such systems to be positive, we require a 
konus density that falls of as r~ p at large radii, where p > 4. 
Furthermore, since we use the criterion that 



dp 
dv 



<9p g al + ^ *9p k on < Q 



(3.1) 



dv dv 
and since for the perfect oblate spheroid 

lim tt-(\, v) oc \- 2 =r~ 4 (3.2), 

A^oo dv 

we require a konus density whose derivative with respect to 
the prolate spheroidal coordinate v decays faster than r~ 4 . 

The specific konus density introduced by GB (their 
equation 19) does not satisfy this criterion. We therefore seek 
another konus density. Kochanek & Rybicki (1996; here- 
after KR) have developed a scheme to derive konus densities 
with a large variety of properties through the introduction 
of 'semi-konus densities'. A semi-konus density is a density 
that is non-zero only in one of the two halfs of the cone of 
ignorance; either the cone at positive or at negative z (here 
z denotes the symmetry axis in Fourier space). Semi-konus 
densities are complex functions and KR showed that the 
real part of any semi-konus density is a konus density. Fur- 
thermore, semi-konus densities have the nice property that 
they are closed under ordinary multiplication: the product 
of two semi-konus densities is itself a semi-konus density. In 
particular any power 



p n (R,z) =K 



Psk(-R, 



(3.3) 



where 5ft denotes the real part, is a konus density as long 
as p sk (R,z) is a semi-konus density. We will refer to konus 
densities of the form (3.3) as konus densities of order n. 

We have experimented with several konus densities and 
finally came up with, what we will refer to as, the general- 
ized konus density (see Appendix B). The generalized konus 
density is characterized by three parameters: a scalelength 
b, an order n, and an additional free parameter /i. Increas- 
ing the latter two, increases the number of oscillations, as 
well as the power p of the decline at large radii. In order 
for the decline of both the density and the derivative (3.1) 
to be sufficiently fast we need to go to order n > 6. This 
is independent of jj,. Throughout this paper we will mainly 
focus on the generalized konus densities with p, = 1. A con- 
tour plot of this konus density with n = 6 in the meridional 
plane is shown in Figure 2. 

4 4 UNCERTAINTIES IN THE CENTRAL 
DENSITY 

We now investigate the maximum konus density one can 
add to or subtract from a perfect oblate spheroid under the 
criterion of equation (2.1). We therefore seek the value 
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Figure 2. Contours of the 6 th order generalized konus density 
p k (R, z 1 45°) for fi = 1. This density distribution projects to zero 
surface brightness for i < 45°. Positive contours are solid, nega- 
tive contours arc dashed. 



where 



f± 



mm 



dv 



J dv 







(4.2) 



Here p£ is the n th order of the generalized konus density (see 
Appendix B). So /+ describes the maximum konus density 
one can add, and /_ the maximum that can be subtracted. 
They are the minimal values of / for which the derivative of 
the total density with respect to v equals zero. This value 
has to be searched over the entire (A, z^)-space. Once /max is 
found we calculate 

/max p k (0,0|i) 



A„ = 



(4.3) 



/max = max(/+,/_), 



(4.1) 



Pgal(0,0) 

The parameter A„ therefore expresses the ratio of the cen- 
tral density of the maximum konus density over the central 
density of the perfect oblate spheroid. 

The results for the generalized konus density of order 
6 with /i = 1 are presented in Figure 3, where we plot the 
logarithm of A„ as a function of the ratio b/a. Here b is 
the scalelength of the konus density and a = \/—a is the 
break radius of the perfect oblate spheroid (see Appendix 
A). The results are shown for four different inclination angles 
(i = 30°, 45°, 60°, and 80°), and two different flattenings q 
(q — 0.8 and 0.5.) For small inclination angles the central 
density of the maximum konus density is of the same order of 
magnitude as the central density of the underlying spheroid 
(i.e., A„ w 1). However, that is only true if the scalelength 
of the konus density is approximately equal to the break 
radius of the perfect oblate spheroid (i.e., b/a ~ 1). For 
higher inclination angles, one can add less konus density to 
the model, reaching / max = for i = 90° . We found that for 
i = 80° , the central density of the maximum konus density is 
no more than ~ 10% of p ga i(0, 0), and only if the scalelength 
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Figure 3. The logarithm of the parameter A„ as a function of the ratio b/a for the 6 th order generalized konus density (p = 1). The 
results are presented for four different inclination angles: i = 30° (solid lines), i = 45° (dotted lines), i = 60° (short dashed lines), and 
i = 80° (long dashed lines). The panel on the left is for a perfect oblate spheroid with an ellipticity e of 0.2, whereas the results in the 
panel on the right are for a model with e = 0.5. 
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Figure 4. The same as Figure 3, except that we now plot the 
maximum konus densities for different orders n of the generalized 
konus density and for different values of fi. All lines are for an 
oblate spheroid with e = 0.2, seen under an inclination angle of 
60°. The thick, solid line is for /i = 1 and n = 6. The dashed 
lines are for /i = 2 and /i = 4, and the dotted lines are for higher 
order; n = 7 and n = 8. As can be seen, central density of the 
maximum konus density increases with both fi and n for large 
scalelengths b. However, for small scalelengths (b/a < 1), A v is 
highest for /i = 1, n = 6. 



of the konus density becomes small. In Figure 4 we show how 
increasing /i or going to higher order n affects the maximum 
konus density. As can be seen, at b/a > 1, the central density 
of the maximum konus density increases with both /i and n. 
For small scalelengths b, however, A„ is maximal for the 



konus density with smallest /x and n (i.e., /x = 1, n = 6). In 
fact the absolute maximum of A„ is reached for the konus 
density with the smallest allowed values of fi and n. 

4.1 4.1 The luminosity profiles 

Any konus density Pk(R, z\io) projects to zero surface bright- 
ness when seen under an inclination angle i < to. However, 
when a galaxy to which we add such a konus density is seen 
more edge-on (i > io), it does not remain invisible. In order 
to examine how realistic galaxy models are with a maximum 
konus density, we project a perfect oblate spheroid with a 
maximum konus density Pk(R, z|30°) under an inclination 
angle of 90° (i.e., we observe this galaxy edge-on), and com- 
pare the luminosity profiles along the major and minor axis 
with and without the maximum konus density. The results 
are shown in Figure 5, where we plot the logarithm of the 
projected density versus the logarithm of the radius from 
the center projected on the sky R. We plot the luminosity 
profiles along the major and the minor axis. We consider the 
/i = 1, n = 6 konus density with b/a = 1.3. For this ratio of 
the characteristic lengths a maximum value of A„ was found 
(see Figure 3). 

As can be seen the presence of the (maximum) konus 
density is clearly visible. Its presence introduces wiggles in 
the luminosity profile along the minor axis. Galaxies with 
wiggles in their luminosity profiles are known. However, in 
all such cases the wiggles are mainly visible along the major 
axis profiles. Examples of such cases are galaxies that have 
embedded nuclear disks (see e.g., Ferrarese et al. 1994, van 
den Bosch & Jaffe 1996). The konus densities used here how- 
ever, reveal themselves by the presence of wiggles along the 
minor axis only. Since no such galaxies are known to exist, 
we can use this as an empirical criterion on the maximum 
konus density. We can make this criterion more qualitative 
by defining a maximum amplitude of the wiggle along the 
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Figure 5. The luminosity profiles along the major and the minor axis of a perfect oblate spheroid with intrinsic flattening q = 0.8 seen 
edge-on (dashed lines). In addition, we also plot the luminosity profile of the same spheroid but to which we have added the maximum 
konus density p^(R, z|30°). The scalelcngth of the konus density b = 1.3a, where a is the break radius of the spheroid. We have plotted 
the logarithm of the projected surface brightness (in arbitrary units) vs. the logarithm of the radius on the sky R in units of a. 




Figure 6. The maximum amount of konus density K v that can be added to a perfect oblate spheroid with q = 0.8, as determined from 
three different criteria: p > (solid lines), dp/dv < (dashed lines), and the empirical criterion that the minor axis surface brightness 
profile should not have wiggles (dotted lines). The results are shown as function of b/a for two different inclination angles (i = 30°, and 
60°). 



minor axis luminosity profile. We define the maximum konus 
density according to the empirical criterion, as that konus 
density for which 

max |^ ga i - jiitotl = 0.05. (4.4) 

Tt. 

i.e., for which the maximum difference between the lumi- 
nosity profiles with (/x to t) and without (/i ga i) konus density 
is 0.05 magnitudes. Given the photometric accuracy with 
which luminosity profiles can be measured observationally, 



any wiggle with an amplitude bigger than 0.05 magnitudes 
should be detectable. In Figure 6 we plot the maximum 
konus density A„ that can be added to a perfect oblate 
spheroid with q — 0.8 as derived by using three different cri- 
teria: p > (solid lines), dp/dv < (dashed lines), and the 
empirical criterion (4.4) that the minor axis surface bright- 
ness profile, when projected edge-on, should not have wiggles 
bigger than 0.05 magnitude (dotted lines). The results are 
shown as a function of b/a for two different inclination angles 
(i = 30°, and 60°). In almost all cases the dp/dv < crite- 
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Figure 7. The ratio A„/Ag as function of the ratio b/a of the 6 th order generalized konus density (p = 1). The results are presented for 
three different inclination angles: i = 30° (solid lines), i = 60° (short dashed lines), and i = 80° (long dashed lines). The panel on the 
left is for a perfect oblate spheroid with an ellipticity of 0.2, whereas the results in the panel on the right are for a model with e = 0.5. 



rion is dominant. Only for low inclinations, and b/a ~ 1 — 3 
is the empirical criterion more strict. 



5 5 AN APPROXIMATE CRITERION 

Observations with the Hubble Space Telescope (HST) have 
revealed that all elliptical galaxies have central cusps in their 
density distribution (e.g., Ferrarese et al. 1994; Lauer et 
al. 1995). The perfect oblate spheroid we have studied so 
far has a density distribution that becomes constant inside 
its break-radius a = \J —a, and therefore is not a realistic 
model for an elliptical galaxy. The advantage, however, of 
the perfect oblate spheroid is that it is a Stackel potential 
so that we can use the criterion dpjdv < 0. Unfortunately, 
no cusped Stackel models are known. In order to investigate 
the uncertainty in the central density of realistic, centrally 
cusped, elliptical galaxies, we therefore need to use another 
criterion. 

We use the criterion that the derivative of the total 
density distribution along a circular curve with radius r = 
\/ R 2 + z 2 has to be negative: 



where 6 = arctan(2/7?). We examine the accuracy of this ap- 
proximate criterion, by comparing the values of A that we 
find for the perfect oblate spheroid using both the proper 
criterion (A„) and the approximate criterion (Ag). The re- 
sults are shown in Figure 7, where we plot the ratio A„/A# 
as a function of b/a for two different Battenings of the under- 
lying perfect oblate spheroid (g = 0.8, and 0.5), and three 
different inclination angles {i — 30°, i = 60°, and i = 80°). 
The approximate criterion always leads to an overestimate of 
the maximum konus density. The error is larger for smaller 
scalelengths b of the konus density, for smaller inclination 
angles, and for less flattened galaxies. 



Having shown that the criterion (5.1) leads to reason- 
ably accurate approximations of the maximum konus den- 
sity for the perfect oblate spheroid, by no means ensures that 
it remains reasonable for, for instance, cusped galaxies. In 
fact, the finding that the approximate criterion leads to very 
accurate maximum konus densities when the scalelength is 
large, is somewhat trivial given that at large radii the prolate 
spheroidal coordinate v becomes equal to 8. Nevertheless, 
since the approximate criterion enforces a certain smooth- 
ness on the density distribution, we feel confident that this 
criterion will lead to more accurate maximum konus densi- 
ties than simply using the criterion that the density has to be 
positive. We showed in Section 2, that under that naive cri- 
terion, one will overestimate the maximum konus densities 
(see also Fig. 6). In Section 7 we show that the approximate 
criterion leads to similar maximum konus densities as the 
empirical criterion (4.4). 

6 6 CUSPED DENSITY DISTRIBUTIONS 

We now investigate the maximum konus density one can 
add to a cusped density distribution. We consider the so 
called (a, /3)-models, which have proven to be useful repre- 
sentations of elliptical galaxies (e.g., Qian et al. 1995). The 
density distribution, which is given by 

= ( 1+ (^)) ' W 

is stratified on concentric ellipsoids m = ^/ R? + z 2 /q 2 of 
constant flattening q. For R 3> R b they fall off as m°* +2/3 , 
while for R <JC Rt they have p oc m a . We will restrict our- 
selves to models that have a + 2/3 = —4 so that at large radii 
their projected surface brightness falls off as R~ 3 . 

We use the approximate criterion to investigate the 
dependence of the maximum konus density we can add 
to the density distribution of equation (6.1), as function 
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Figure 8. The ratio of the central density of the maximum konus density over the density of the (a, /3)-model at its core-radius Rt as 
function of the ratio b/Rb- Results are plotted for four different values of the cusp steepness a, and for three different inclination angles: 
i = 30° (solid lines), i = 60° (dashed lines), and i = 80° (dotted lines). The maximum konus density increases strongly with increasing 
cusp steepness and decreasing konus-scalelength 6. 



of the cusp steepness a. We consider a strongly flattened 
galaxy (q = 0.5) seen under three different inclination an- 
gles (i — 30°, i — 60°, and i = 80°). When the cusp steep- 
ness a < the central density of the (a, /3)-model is infinite. 
We can therefore not calculate the uncertainty of the cen- 
tral density Ae. Instead, we calculate the ratio A' e of the 
central density of the maximum konus density over the den- 
sity po of the (a, /3)-model. We restrict ourselves to a konus 
density with fi = 1 and n = 6. The results are shown in 
Figure 8. When no cusp is present (i.e., a = 0) we again 
find a maximum A' e at b/Rb ~ 1. However, for density dis- 
tributions that are cusped we find that A' g keeps increasing 
when b/Ri, — > 0. In the following section we will investigate 
the importance of this for galaxies where strong indications 
are found for the presence of a nuclear BH by discussing a 
specific example; the compact elliptical M32. 



7 7 APPLICATION TO M32 

To illustrate the importance of konus densities to real galax- 
ies we here apply our analysis to M32. Strong evidence ex- 
ists that M32 harbors a nuclear BH of ~ 1.8 x 10 6 M (van 
der Marel et al. 1994, Qian et al. 1995, Dehnen 1995). This 
evidence is based on the observations of a strong central 
increase of the velocity dispersion. In this Section we in- 
vestigate whether redistributing the central mass of M32 
according to the konus density distribution can account for 
this increase without adopting the presence of a nuclear BH. 

Although the inner part of the galaxy can accurately 
be described by an (a, 13) -model, this does not fit the outer 
parts of the galaxy. Multiplication of the (a,/3)-model with 
an additional factor can adequately describe the strong de- 
crease of the density beyond ~ 20" . We therefore follow the 
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Figure 9. The logarithm of the parameter as a function 
of the scalclcngth b of the generalized konus density for the 
(a, /3,7)-model of M32. The lower two curves are for the konus 
density with fi = 1 and n = 6, whereas the upper two curves cor- 
respond to a fi = 4, n = 6 konus density. Solid curves correspond 
to criterion (5.1), whereas the two dashed curves show the maxi- 
mum konus density as derived from the empirical criterion (4.4). 
Both criteria give values for A of the same order of magnitude, 
with the empirical criterion being slightly more strict. 

approach taken by van der Marel et al. (1996), and describe 
the density distribution of M32 by the following formula 
(hereafter referred to as a (a,/3,7)-model) 



Major axis 
1 rn — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — r 
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Projection of this density distribution and fitting to the 
observed surface brightness (from both HST and ground- 
based data) results in a good fit to the data for the parame- 
ters: a = -1.435, fi = -0.423, 7 = -1.298, R b = 0.55", 
and R c = 102.0" (van der Marel, priv. communication). 
Throughout we will assume a distance of M32 of 700 kpc, 
so that 1" corresponds to 3.4 pc. 

An important problem with M32 is the fact that the 
inclination angle is unknown. Van der Marel et al. (1994) 
constructed axisymmetric models of M32 for different as- 
sumptions of the inclination angle (i = 90° and i = 50°). 
After solving the Jeans equations they concluded that both 
models require the presence of a ~ 1.8 x 10 6 Mq BH in or- 
der to fit the observed velocity dispersions in the center. If 
indeed M32 is observed edge-on there is no uncertainty in 
the deprojection of the surface brightness so that no konus 
density can be added. Here we will focus on the i = 50° 
model for which van der Marel et al. (1994) derived a mass- 
to-light ratio in the Johnson V-band of 2.36. Given that the 
observed flattening of M32 is equal to 0.73, this results in an 
intrinsic flattening of the ellipsoids m of 0.452 (if i = 50°). 
The density p of equation (7.1) is 1.76 x 10 5 M Q pc" 3 . 

We have calculated the maximum konus density p£ that 
we can add to this (a,p,j)-mode\ of M32. Since at large 
radii m oc m a+2/3+27 ~ m~ 4 ' 8 we again need n > 6. We 
use the approximate criterion dp/d9 < 0, and calculate the 
parameter A' g for different values fi, n, and b of the konus 
density. The results are shown in Figure 9 for two different 
values of n (fj, = 1 and fi = 4 both with n = 6). Unlike with 
the perfect spheroid, A is larger for higher (iorn for each 
value of the scalelength b (compare Figure 4). In principle, 
we can make A' arbitrarily large by going to small b and 
large n or p. 

In addition, we have determined the maximum konus 
density allowed according to the empirical criterion (4.4) 
that the edge-on projected surface brightness along the mi- 



Minor axis 




Log(R) [arcsec] 



Figure 10. The logarithm of the edge-on projected surface brightness (in arbitrary units) of the (a,/3,7)-model of M32, both with (solid 
lines) and without (dashed lines) the maximum konus density (/j, = 1, n = 6) derived from the approximate criterion. 
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Figure 11. The maximum mass fraction r max (equation 7.4) as 
function of /i + n. Since /j. > 1 and n > 6 we have a minimum 
/i+n = 7. The solid line shows I\ on for increasing order n (fi = 1), 
whereas the influence of increasing fj, (with n = 6) is indicated by 
the dashed line. 

nor axis should not have wiggles exceeding 0.05 magnitude 
(dashed lines in Figure 9). Both criteria give rather similar 
values for A, although the empirical criterion is more strict. 
This is evident from Figure 10, where we plot the luminosity 
profiles for the (a,/3,7)-model of M32 both with (solid lines) 
and without (dashed lines) the maximum konus as derived 
from the approximate criterion. 



7.1 7.1 Dynamical significance of konus densities 

We have shown that by increasing /i or n we can make A' 
arbitrarily large. However, A' is the central konus density 
divided by the density of the galaxy at m = R^. Since the 
central density of a cusped galaxy is by definition infinite, 
it is not a priori clear what the real, dynamical significance 
of these maximum konus densities is. In order to compare 
the dynamical significance of the (maximum) konus density 
with that of the inferred BH, we have calculated the ratios 



(r) 



Mg a i(r) +M BH 
M g ai(r) : 



and 



■ koj 



M ga i(r) + M kon (r) 
M g ai(r) 



(7.2) 



(7.3) 



Here M(r) is the mass inside a sphere of radius r, and M ga i 
is the mass of the (a,/3,7)-model. Although the total mass of 
the konus density is zero, at sufficiently small radii there is 
only a positive density contribution resulting in a non-zero 
positive Mkon(r). Since the central konus density is finite, 
TkonM is zero at r = 0, reaches a maximum at a certain 
intermediate radius, and decreases to zero again at larger 
radii. This (r) depends on the scalelength b of the 

konus density, as well as on p and n. For each p and n we 
calculate the maximum mass fraction, defined as 



This maximum mass fraction as function of pi + n is shown 
in Figure 11. Although increasing p + n increases the central 
density of the maximum konus density (see Figure 9), the 
number of oscillations increases as well. Therefore, the cen- 
tral region where pk > becomes smaller and we find that 
r m ax decreases with increasing p + n. For our generalized 
konus density we find an absolute maximum mass fraction 
of ~ 1.04 only. This maximum is reached for a konus density 
with b = 0.22", p = 1, and n = 6. It is not only important 
what the absolute mass fraction is, but also at what radius 
it is reached. In Figure 12 we have plotted the radial run 
of Tbh- If r max = 1.04 is reached at 5 pc, it is compara- 
ble to Fbh- However, the maximum mass fraction occurs at 
very small radii (r as 0.7pc), where Fbh is many orders of 
magnitude larger. If the scalelength b of the konus density 
increases, the radius where Fkon(r) is maximal increases, but 
since A decreases strongly with scalelength (Figure 9), the 
maximum Fk on also decreases with b. We therefore conclude 
that the dynamical influence of the maximum generalized 
konus density is negligible compared to that of the inferred 
BH. 

7.2 7.2 Cusped konus densities 

The Fourier Transform of the generalized konus density de- 
cays exponentially at large frequency k (see Appendix B). 
This enforces dpk/dr to go to zero at small radii r. Here 
we investigate whether konus densities with a central den- 
sity cusp can add more mass to the center. In order for a 
konus density to be cusped, i.e., pk oc r~ a (0 < a < 3) at 
small radii, its Fourier Transform should decay as fc a ~ a at 
large frequencies k. The requirement a < 3 ensures that the 
density cusp has finite mass. We can constrain this regime 
even further by taking into account the fact that every konus 
density projects to zero surface brightness when projected 
face-on (i.e., with i = 0°). Since 



Jo 



Mr = co, 



(7.5) 



Tmax — max 



maxTkonfr) 



(7.4) 



for a > 1 and a > 0, the projected surface brightness at the 
center can only be equal to zero if a < 1. Therefore, konus 
densities can at most be moderately cusped, with a central 
density gradient less steep that r _1 . 

As we have seen in the previous section, maximizing 
Tkon, not only requires high central density, but also a suf- 
ficiently large area in the center where pk > 0. A cusped 
konus density with this property, in order to have zero to- 
tal mass, must have regions outside the center with strongly 
negative density. Addition of such a konus density, is there- 
fore likely to result in luminosity profiles (when projected 
edge-on) with large wiggles. In other words, it seems likely 
that the empirical criterion (4.4) becomes most strict. In 
fact, already for our generalized konus density, this criterion 
is the most restricting (see Figure 9). In order to quantify 
this, we have tried to construct cusped konus densities, but 
were unable to find an analytical pk(-R, z) that obeys the 
physical restrictions of konus densities, i.e., zero total mass, 
sufficiently rapid decay at large radii, and no discontinu- 
ities. We therefore constructed a density distribution that 
resembles a cusped konus density. Consider the density dis- 
tribution 

pc{R,z\i) = 
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Figure 12. The total mass inside a sphere of radius R, Mtot(R), 
over the mass of the M32 density model inside the same sphere 
as function of radius R. The solid line shows the case in which 
the total mass is build up of the (a, /3,7)-modcl of M32 with in 
addition a BH of 1.8 X 10 6 Mq. The dashed line is for the case 
in which we have added the maximum generalized konus density 
with b = 1.0" to the density model of M32. 
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Figure 13. A contour plot of the cusped 'konus '-density 
pc (R, z 1 50° ) for a = 0.25. The parameters b and d are taken 
such that p c (R,z\50°) obeys the criteria (7.7) and (7.8). 



p Q r- a e- r/c e-^tanh 



WR 2 +d 2 tan" 1 ' 



-)■ ( 7 - 6 ) 



Here r = V-R 2 + z 2 , and a, fe, c, and d are free parame- 
ters. This density is cusped in the center (with power law 
slope a), and decays exponentially at large radii r » c. At 
z — \/R 2 + d 2 tan -1 i the density changes sign, becoming 



negative at larger z. We adopt a and c as free parameters 
and use the physical properties of genuine konus densities 
to determine b and d. Since we have added the factor e~ R / b , 
which equals unity for z — 0, we can use the requirement 
that the line-of-sight integral along the minor axis should 
yield zero, i.e., 



p c (R = 0, z\i)dz = 0, 



(7.7) 



to numerically evaluate the ratio d/c. Once d is known, we 
determine the ratio b/c, by numerical evaluation of the root 
of 

poo poo 

= 4tt / dRR dz p c (R,z\i) = 0. (7.8) 
Jo Jo 

A contour plot of a cusped 'konus'-density with a = 0.25 
that obeys these criteria is shown in Figure 13. Although 
(7.6) is not a konus density, we can interpret it as such and 
use the empirical criterion (4.4) to determine the maximum 
'konus'-density. The empirical criterion uses the edge-on pro- 
jected surface brightness along the minor axis, and along this 
axis (7.6) at least resembles a genuine cusped konus density, 
in that it obeys criterion (7.7). As we have seen for the gener- 
alized konus density, the maximum mass fraction decreases 
for an increasing number of density oscillations. Our cusped 
'konus'-density has a minimal number of density oscillations, 
and is therefore biased towards larger Fkon- We have calcu- 
lated, as function of a and scalelength c, the max r F^ OI1 (r ) 
for our cusped 'konus'-density when added to M32. The re- 
sults are shown in Figure 14. Similar as for the generalized 
konus density, (r) increases with decreasing scale- 

length. More importantly, increasing the cusp slope a de- 
creases Fkon- This is due to the fact that for larger a the 
konus-criterion (7.7) ensures that d decreases, so that the 
central region where p c > becomes smaller. Evidently, this 
effect is stronger than the increase in cusp slope. We there- 
fore conclude, that even cusped konus densities have a negli- 
gible dynamical effect on the central region of M32, as com- 
pared to that of the inferred BH mass. In the analysis above 
we have used a BH mass of 1.8 x 10 6 M Q . More recent studies 
hint towards a more massive BH of ~ 3 x 10 6 M (Bender, 
Kormendy & Dehnen 1996; van der Marel et al. 1996). This 
only strengthens the conclusions outlined above. 



8 8 CONCLUSIONS 

The deprojection of axisymmetric systems is indeterminate 
when the inclination angle i ^ 90° . To any such system one 
can add densities which, after projection, result in zero sur- 
face brightness. Since these density distributions have zero 
total mass, addition or subtraction of konus densities is sim- 
ilar to redistributing the mass of the system. In this paper 
we have investigated to what extent one can redistribute the 
mass in elliptical galaxies. In particular we investigated the 
uncertainties of the central densities of such galaxies due to 
the non- uniqueness of the deprojection. 

Since konus densities have regions with both positive 
and negative density, there is a limit on how much konus 
density can be added to any density distribution while main- 
taining positivity. Although this criterion will result in a 
maximum konus density, it ignores whether a model which 
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Figure 14. The maximum of I\ on over radius r for our cusped 
'konus'-density p c (7.6) as function of its scalclength c, and cusp 
slope a. Increasing the cusp slope, decreases the mass fraction 
added to the center by the 'konus'-density, therewith resulting in 
a smaller dynamical effect as compared to that of the BH. 

such a maximum konus density can actually be built from 
its orbit building blocks. We therefore use a more realistic 
(and more strict) criterion to determine the maximum konus 
density. 

We used both a physical and an empirical criterion. The 
physical criterion is based on the principle that the density 
should be monotonically decreasing along curves that fol- 
low thin tube orbits. This criterion, although sufficient is 
not strict. For Stackel potentials, these thin tube orbits are 
described by the curves of constant prolate spheroidal co- 
ordinate A, and the criterion can be made quantitative. For 
any other potential, no analytic description of the curves 
is known. In those cases we use an approximate criterion 
that states that dp/ 88 < along circles of constant radius 
r = \/ R 2 + z 2 with 8 = arctan (z/i?). The empirical cri- 
terion we used is based on the fact that adding too much 
konus density will result in wiggles along the minor axis lu- 
minosity profile, when the galaxy is seen edge-on. Since no 
such systems have ever been observed this puts an empirical 
constraint on the maximum konus densities. 

We calculated the ratio of the central density of the 
maximum konus density to the central density of the perfect 
oblate spheroid to which that konus density can be added. 
For sufficiently small inclination angles one can add konus 
densities whose central density is comparable to that of the 
perfect oblate spheroid. For highly inclined galaxies the un- 
certainty of the central density is small however: < 10% for 
i 80° . Perfect oblate spheroids have a break radius inside 
which the density is constant. Elliptical galaxies, however, 
are known to harbor central density cusps. We therefore 
applied the approximate criterion to cusped, axisymmetric 
models that have p oc r~ a at small radii. The steeper the 
cusp, the more konus density one can add. The central den- 
sity of the maximum konus density increases strongly with 
decreasing scalelength of the konus density and with decreas- 
ing inclination angle. 



In order to better understand the dynamical signifi- 
cance of maximum konus densities, we applied our analy- 
sis to the specific example of M32; a compact elliptical for 
which strong indications are found for the presence of a cen- 
tral BH of ~ 2 x 10 6 M Q . The central density distribution 
of M32 is strongly cusped (p oc r _1 ), and we used both 
the approximate criterion and the empirical criterion to de- 
termine the maximum amount of konus density. The latter 
was found to be more strict. Although for sufficiently small 
scalelengths of the konus density its central density can be 
very large, its dynamical influence on the central region of 
M32 is negligible. This is due to the cusped nature of the 
central density distribution of M32. 

We have also investigated the dynamical influence of 
cusped konus densities. We have shown, based on a simple ar- 
gument, that konus densities can only be moderately cusped 
(pkon(f) oc r~ a with a < 1). We were unable to find an an- 
alytic, cusped konus density, and therefore used a cusped 
density distribution that resembles a cusped konus density. 
We made sure that this density distribution has zero total 
mass, and when projected face-on, has zero surface bright- 
ness in the center. Using the empirical criterion, we showed 
that for increasing cusp slope a, the mass fraction that the 
konus adds to the center decreases. This is due to the fact 
that an increase in cusp slope, results in a decrease of the 
central region where pkon > 0. 

We conclude therefore that konus densities play only a 
marginally important role in the study of the dynamics of 
the central regions of (cusped) elliptical galaxies. Although 
considerable amounts of konus density can be added, when 
its scalelength is sufficiently small, its dynamical influence 
on the central dynamics remains negligible. Especially, we 
conclude that the uncertainties in the central mass of ax- 
isymmetric galaxies due to the non-uniqueness of the de- 
projection are sufficiently small that they cannot exclude 
the requirement of a central BH based on the finding of too 
small central densities in these galaxies. This is more and 
more true in galaxies seen under larger inclination angles. 

A major shortcoming in the investigations discussed in 
this paper is the fact that only one particular class of konus 
densities was used. However, the different free parameters of 
this generalized konus density allowed a wide range of konus 
densities to be investigated. Although the cusped 'konus' 
density distribution was not a real konus density, it at least 
obeyed two konus criteria, and we merely used it to inves- 
tigate the changes in mass fraction when cusp slope is in- 
creased, rather than to interpret the absolute mass fraction 
for this 'konus' density. We are therefore confident that the 
conclusions reached in this paper are at most very moder- 
ately dependent on our actual choice of konus densities used. 
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APPENDIX A: THE PERFECT OBLATE 
SPHEROID 

Kuzmin (1953; 1956) showed that their exists exactly one 
oblate potential for which the density is stratified on similar, 
concentric aligned spheroids, and for which the equations of 
motion are separable. This model has become known as the 
perfect oblate spheroid (de Zeeuw & Lynden-Bell 1985), and 
has a potential of Stackel form, i.e., 



(A + 7)G(A)-(i/ + 7 )G(i/) 
X-v 



(Al) 



where G(t) is an arbitrary function (r = A, v), and (A, (f>, v) 
are prolate spheroidal coordinates, in which the equations of 
motion are separable. The coordinates A and v are denned 
as the roots for r of 

?2 2 



R 1 



0, 



(A2) 



r + a t + 7 

where R and z are the cylindrical coordinates, a and 7 are 
negative constants, and we choose —7 < v < —a < A. The 
foci of the coordinates are at (R, z) = (0, ±A), with A 2 = 
7 — a. The coordinates (R, z) can be written as 



D 2 {\ + a)(v + a) 2 (A + 7)^ + 7) 



(A3) 



a — 7 7 — a 

The density distribution of the perfect oblate spheroid 



is given by 
p{R, z) = po 



(A4) 



(a 2 + rn 2 ) 2 ' 

where a = \/—a, and m 2 = R 2 + z 2 /q 2 are the spheroids 
with flattening q = ^/7/a. The density distribution in the 



equatorial plane has a break radius at a. At radii r » a 
the density falls of as r~ 4 (consistent with the fact that the 
surface brightness of most galaxies falls off as approximately 
R~ s ). At small radii (r << a) the density is constant and 
equal to the central density po ■ The density distribution can 
be expressed in prolate spheroidal coordinates as 



p(X, v) = p 



»7 



(A5) 



APPENDIX B: 
DENSITY 



A GENERALIZED KONUS 



Here we follow the approach taken by KR, and consider the 
class of scmi-konus densities which can be written as 



p sk (R, z\i) 



[/.(>+) + /.(>-)] 



.,[/c(a + )-/c(q-)] ;(B1) 



where a± = \/ R 2 + z 2 cosi±z, and f s (x) and f c (x) are odd 
and even functions, respectively, defined by 



fc(x) 



and 



/.(*) 



1 

4^2 



OC 

4^/ 



dfc g(k) cos(kx) 



dk g(k) sin(fcr). 



(B2) 



(B3) 



Upon choosing a weight function g(k) oc fc M 1 e bk one de- 
rives with the help of formulae 3.944.5 and 3.944.6 in Grad- 
shteyn & Ryzhik (1980): 

M x ) = m cos A* arctan (f) , (B4) 



and 



(b 2 + x 2 y/ 2 



p arctan (^) 



(B5) 



where F(a;) is the gamma function. These functions define a 
continuous class of konus densities pk(R, z\i) = 5ftp s k(-R, z\i) 
for p > 0. Here 5ft/ denotes the real part of /. For p = 1 
the semi-konus density (Bl) reduces to the first order semi- 
konus density given by KR in their equation (20). In the 
following we will only consider integer values of p. 

In the limit r — V R 2 + z 2 — » 00 the konus densities 
Pk(R, z\i) decline as r~' M+1 ^ for p odd and as r^' M+2 - ) for p 
even. This continuous set of konus densities therefore can be 
chosen to decline arbitrary rapidly at large radii. In Figure 
15 we show contour plots for the cases with p — 1, 2, 3, and 
4. Note that for p > 3 these konus densities also change sign 
in the equatorial plane. 

As we prove in Appendix C, the generalized konus den- 
sity Pk(r, z\io) with p — 1 has infinite total mass. For i < io 
it projects to constant, but non-zero surface brightness. 
Therefore, addition of such a konus density is no longer 
equivalent to simply redistributing the density. For p = 2 
the total mass is finite but positive, whereas for p > 2 the 
total mass is zero (see Appendix C). 

As was shown by KR, any function 



p n (R,z) = 5ft 



plk(R, 



(B7) 
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Figure 15. Contours of the generalized pk(R, z\60°) for p, = 1,2,3, and 4. Each of these density distributions projects to zero surface 
brightness for i < 60°. Positive contours are solid, negative contours are dashed. 



is a konus density as long as pl k (R, z) is a semi-konus den- 
sity. KR have presented examples of such higher order konus 
densities for the case of p, — 1 of our generalized konus den- 
sity. 

As discussed in the section 3, in order to ensure that re- 
alistic axisymmetric galaxy models can be constructed with 
the addition or subtraction of konus densities, one needs 
konus densities that decline sufficiently fast. Any konus den- 
sity of the form (Bl), with f s (x) odd (f' 3 (0) / 0) and 
f c (x) even (/ c (0) 7^ 0) has pj? oc r~ n for n even, and 
dp%:/d9 oc r~' n ~^ for n odd along lines a± = 0. Here 
9 = arctan(.z/.R), but the same is true for the derivative 
with respect to the prolate spheroidal coordinate v. In order 
to have a konus density that declines faster than r~ 4 at large 



radii, and whose derivative dp k /d9 declines faster than r~ 4 , 
we need to go to order n > 6. Besides a more rapid asymp- 
totic decline for higher order n, the konus density will also 
show a greater number of sign alterations. 



APPENDIX C: THE TOTAL MASS OF THE 
GENERALIZED KONUS DENSITIES 

Let H(k) be the Fourier Transform of the semi-konus den- 
sity p s k(r, z\9), where 9 — 90° — i. We will use spherical 
coordinates (fc, tp, <f>) for the Fourier space vector k. Then 
H(k) — H(k,tp) since the Fourier Transform of the semi- 
konus density is symmetric around the k z axis. We follow 
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Figure 16. The build-up of the Fourier Transform H(k, tp) of the 
scmi-konus density inside a cone with half-angle 8, by a weighted 
sum of spheres. At each point (k, tp) with tp < 8 there are two 
spheres with radii fc+ and k- that contribute to H{k,tp). 



the approach taken by KR and build up H(k) by a weighted 
sum of kernel functions f(k, ko), i.e., 



oo 

H(k) = J 



dk g(k )f(k, ko). 



For our specific semi-konus density of the form (Bl), the 
kernel is a shell of radius | ko | = ko shifted along the k z axis 
by an amount &o/sin#, and 



f(k,ko) = 



4 


k — ko/ 


sir 







ko 


] 




ko 





(C2) 



where 8(x) is the delta function. This is illustrated in Figure 
16, where we have plotted the cone with half-angle in the 
meridional plane of k. As can be seen, at any point (k, tp, cp) 
(with tp < 0) there are two shells that contribute weight to 
that point. We define the function f(ko) as 



f(ko) = \ k — fco/sin^ — \ko\ 



k 2 + 



k> 2 - - 2fc-^ cost/; - k , 
sm 2 sin 



(C3) 



From the normalization relation of 5-functions we have that 
5(k - k') 



5[f(ko)] = 



f'(k') 



(C4) 



where k' are those fco for which /(fco) = 0. Solving /(fco) = 
we find two solutions k' = k± 



k± = k 



sin 
cos 2 *? 



cos tp ± ^/ cos 2 l/> — COS 2 



(C5) 



Calculating the derivative /'(fco) and using equations (C2), 
(C4), and (C5) we find (for tp < 0) 



H(k) = 



sin 



\J cos 2 tp — cos 2 



g(k+) g(k-) 



(C6) 



The function H(k) can also be written as the Fourier 
Transform of the semi-konus density, i.e., 



H(k)= J d 3 fp sk (r\0) exp(-ifc-r). 

Therefore, one has that 

#(0,0,0) = J d :i r Psk (r\0) = M kon , 



(C7) 



(C8) 



where Mkon is the total mass of the semi-konus density. 

The weight function we used to construct the general- 
ized konus density is g(k) = k^~ 1 e~ bk . Since k± oc k wc 
derive from equations (C6) and (C8) that 



: C ■ lira fc M ~ 2 



-6/(0)* 



-6g(0)fc 



(C9) 



where C is a constant that depends on jj,, b, and 0, and / 
and g are functions of the angle tp. Therefore, for /i = 1 the 
total mass of the konus density is infinite. For fi = 2 one has 
Mkon > 0, but finite, and for /i > 2 the total mass is equal 
to zero. 

We now show that the higher order konus densities 
Pk(r, z\0) (n > 2) of the generalized konus densities have 
zero total mass. Therefore we write the Fourier Transform of 
the semi-konus density p" k (r, z\0) as H£(k), where n denotes 
the order of the semi-konus density, and b the scalelength. 
From equation (C6) it is straightforward to show that 



( C1 ) Hl{ak) = a»- 2 Hl b {k), 



(C10) 



where a > is a constant, and ak — (ak,tp,(p). The FT of 
the second order semi-konus density is the convolution of 
H b (k) with itself, i.e., 



Ht(k) = J d A k' Hl(k-k') Ht(k'). 
Then 

Hl{ak) = J d 3 k' Hl{ak-k') Hl(k'). 



(Cll) 



(C12) 



Since k' is a dummy variable of integration, we can substi- 
tute k' — az to find 



(C13) 



Hi(ak) = a 3 J d 3 zHl(a(k-z)) Hl(az) 
= a 2 ^ 1 J d 3 zHl b {k-z) Hl b {z) 



2/i-l 



H 2 ab (k). 



From equation (C13) we see that whereas Hl(k) oc fc M 1 at 
small radii, we find that H b (k) oc fc 2M_1 at small radii. From 
this analysis it is clear that the total mass M to t = H b (0, 0, 0) 
is equal to zero for n > 2. 
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